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Abstract 

Since the seminal work of Ambrosetti and Prodi, the study of global folds 
was enriched by geometric concepts and extensions accomodating new exam¬ 
ples. We present the advantages of considering fibers, a construction dating to 
Berger and Podolak’s view of the original theorem. A description of folds in 
terms of properties of fibers gives new perspective to the usual hypotheses in 
the subject. The text is intended as a guide, outlining arguments and stating 
results which will be detailed elsewhere. 
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1 Introduction 

When we teach the first courses in calculus and complex or real analysis, a great 
emphasis is given to geometric issues: we plot graphs, enumerate conformal map¬ 
pings among special regions, identify homeomorphisms. Alas, this is far from being 
enough: mappings become too complicated soon. Still, the geometric approach, 
especially combined with numerical arguments, is very fruitful in some nonlinear 
contexts. 

It is rather surprising that some infinite dimensional maps can be studied in a 
similar fashion — one may even think about their graphs! The examples which are 
amenable to such approach are very few, and they elicit the same sense of wonder 
that (the equally rare) completely integrable systems do: one is left with a feeling 
of deep understanding. This text is dedicated to some such examples. 

The interested reader could hardly do better than going through the review 
papers by Church and Timourian (ca, m), which cover extremely well the ma¬ 
terial up to the mid nineties. Their approach is strongly influenced by the original 
Ambrosetti-Prodi view of the problem, which we describe in Section 12.21 In a nut¬ 
shell, the global geometry of a proper function F is studied through certain proper¬ 
ties of its critical set C together with its image F{C), along with the stratification 
of C in terms of singularities. 

This much less ambitious text is mainly an enumeration of techniques and of 
some recent developments, some of which have not been published. We mostly take 
the Berger-Podolak route ([!]) which has been extended by Podolak in and, we 
believe, still allows for improvement. Instead of the critical set, we concentrate on 
the restriction of F to appropriate low dimensional manifolds (one dimensional, in 
the Ambrosetti-Prodi case), the so called fibers. 

Essentially, fibers are appropriate in the presence of finite spectral interaction, 
which roughly states that the function F : X ^ Y splits into a sum of linear and 
nonlinear terms, F = L — N and N deforms L substantially only along a few eigen¬ 
vectors spanning a subspace V C X. The domain splits into orthogonal subspaces, 
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X = H®V and the hypotheses on the nonlinearities are naturally anisotropic. Dif¬ 
ferent requests on H and V yield a global Lyapunov-Schmidt decomposition of F: 
on affine subspaces obtained by translating H, F is a homeomorphism and compli¬ 
cations due to the nonlinear term manifest on fibers, which are graphs of functions 
from V to FI. 

Fibers are also convenient for the verification of properness of F. In particular, 
one may search for folds in nonlinear maps defined on functions with unbounded 
domains, which are natural in physical situations. Fibers also provide the conceptual 
starting point for algorithms that solve a class of partial differential equations, an 
idea originally suggested by Smiley ([31], m) and later implemented for finite 
spectral interaction of the Dirichlet Laplacian on rectangles in [5] . 

An abstract setup in the spirit of the characterization of folds as in or like 
the one we present in Section 01 provide a better understanding of the role of the 
hypotheses in the fundamental example of Ambrosetti and Prodi. Elliptic theory 
seems to be less relevant than one might think, it is just that it provides a context 
in which the required hypotheses are satisfied. 

In Section 2, we present the seminal examples — the Dolph-Hammerstein home- 
omorphisms and the Ambrosetti-Prodi fold — in a manner appropriate for our ar¬ 
guments. Fibers and sheets are defined and constructed in Section 3. A global 
change of coordinates in Section 4 gives rise to adapted coordinates, in which the 
description of critical points is especially simple. A characterization of the critical 
points strictly in terms of spectral properties of the Jacobian DF is given. Also, 
the three natural steps to identify global folds become easy to identify. Further 
study of how to implement each step is the content of Sections 5, 6 and 7. The last 
section is dedicated to some examples. 

The text is written as a guide: we try to convey the merits of a set of techniques, 
without providing details. Complete proofs will be presented elsewhere ([?], [?]). 

Alas, we stop at folds. There are scattered results in which local or global 
cusps were identified: again, the excellent survey [13] covers the material up to 
the mid nineties. So far, the description of cusps seems rather ad hoc. There are 
characterizations m), but they are hard to verify and new ideas are needed. On 
the other hand, checking that maps are not global folds is rather simple, a matter of 
showing for example that some points in the image have more than two preimages. 
A numerical example is exhibited in Section [5.31 


2 The first examples in infinite dimension 

Among the simplest continuous maps between Hilbert spaces are homeomorphisms, 
in particular linear isomorphisms. A second class of examples are folds. 

2.1 Homeomorphisms: Dolph and Hammerstein 

Dolph and Hammerstein 1|15]. [I6j l obtained a simple condition under which non¬ 
linear perturbation of linear isomorphisms are still homeomorphisms. A version of 
their results is the following. 

Start with a real Hilbert space Y and a self-adjoint operator L : X C Y Y 
for a dense subspace X of Y. Let cr{L) be the spectrum of L. 

Theorem 1 Let [—c, c] r\(T{L) = 0 and suppose N : Y ^ Y is a Lipschitz map with 
Lipschitz constant n < c. Equip X with the graph topology, ||a;||jf = ||a;||v + 

Then the map F = L — N:X^Y is a Lipschitz homeomorphism. 
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Indeed, to solve F(x) = y, search for a fixed point of 

Cy-.Y^Y, Cyiz)=N{L-\z))+y 

which is a contraction because the operator L~^ : Y ^ Y has norm less than 1/c 
by standard spectral theory and then the map N o L~^ is Lipschitz with constant 
less than njc < 1 . As usual, the fixed point varies continuously with y. Clearly, F 
is Lipschitz. To show the same for F~^^ keep track of the Banach iteration. 

Notice that the statement allows for differential operators between Sobolev 
spaces. Very little is required from the spectrum of L. Clearly, for symmetric 
bounded operators one should take X = Y. 

2.2 Breaking the barrier: the Ambrosetti-Prodi theorem 

What about more complicated functions? Ambrosetti and Prodi m) obtained an 
exquisite example. After refinements by Micheletti and Manes ([H]), Berger and 
Podolak ([?]) and Berger and Church ([S]), the result may be stated as follows. Let 
n C K" be a connected, open, bounded set with smooth boundary (for nonsmooth 
boundaries, see [33]). Let and be the usual Sobolev spaces and set 

X = F[^{n) n iLo(n) and Y = = L^{Vt). The eigenvalues of the Dirichlet 

Laplacian —A:XcY Y are 

cr(—A) = {0 < Ai < A2 < ... —>■ 00}. 

Denote by (/)i the (L^-normalized, positive) eigenvector associated to Ai and 
split X = Fix © Vx, Y = Hy © Vy in horizontal and vertical orthogonal subspaces, 
where Vx =Vy = the one dimensional (real) vector space spanned by (pi. 

Theorem 2 Let F : X ^ Y be F = L — N, where L = —A, N(u) = f(u), for a 
smooth, strictly convex function / : M ^ R satisfying 

Ran f = {a,b), a < \i <b < X 2 ■ 

Then there are global homeomorphisms C : A —>■ Hy © R and ^ : V —^ Hy © R 
for which F{z,t) = F o (~^{z,t) = {z, —t^). 

Said differently, the following diagram commutes. 

X Y 

Ci i ? 

iLy © R ^ iJy © R 

Functions which admit such dramatic simplification are called global folds. The 
vertical arrows in the diagram above are (global) changes of variables and sometimes 
will be maps, but we will not emphasize this point. 

The original approach by Ambrosetti and Prodi is very geometric ([1|). In a 
nutshell, they show that F" is a proper map whose critical set C (in the standard 
sense of differential geometry, the set of points u G A for which the derivative DF{u) 
is not invertible) is topologically a hyperplane, together with its image F{C). They 
then show that F is proper, its restriction to C is injective and F~^{F{C)) = C. 
Finally, they prove that both connected components of A — C are taken injectively 
to the same component of V — F{C). Their final result is a counting theorem: the 
number of preimages under F can only be 0, I or 2. 

Berger and Podolak ([4]), on the other hand, construct a global Lyapunov- 
Schmidt decomposition for F. For Vx = Vy = (pi), consider affine horizontal 
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(resp. vertical) subspaces of X (resp. Y), i.e., sets of the form Hx + t^i, for a 
fixed t S K (resp. y + Vy, tor y € Hy). Let P : Y ^ Hy be the orthogonal 
projection. The map PFt : Hx —>■ Hy , PFt{w) = PF{w + t(j)i), is a bi-Lipschitz 
homeomorphism, as we shall see below. Thus, the inverse under F of vertical lines 
y + Vy, tor y G Hy are curves : R ^ Vx C X — >• Hx, which we call fibers. Fibers 
stratify the domain X. Thus, to show that F is a global fold, it suffices to verify 
that each restriction F : ay ^ Vy ^ M., essentially a map from R to R, is a fold. 

After such a remarkable example, one is tempted to push forward. This is not 
that simple: if the (generic) nonlinearity / is not convex, there are points in Y with 
four preimages ([?]), so the associated map F : X ^ Y cannot be a global fold (for 
a numerical example, see Section . 

3 Fibers and height functions 

Fibers come up in [4] and [32] for maps associated to second order differential 
operators and in |25j in the context of first order periodic ordinary differential 
equations. Due to the lack of self-adjointness, the construction in [5S| is of a very 
different nature. We follow m and [34] . which handle Lipschitz maps, allowing the 
use of piecewise linear functions in the Ambrosetti-Prodi scenario, namely / given 
by f'{x) = a or b, depending if a: < 0 or a; > 0 ([Hj, [TOji. 

Let X and Y be Hilbert spaces, X densely included in Y. Let L : X CY Y 
be a self-adjoint operator with a simple, isolated, eigenvalue Xp, with eigenvector 
4>p & X with ||(/)p||v = 1. Notice that Xp may be located anywhere in the spectrum 
cr(L) of L. As before, consider horizontal and vertical orthogonal subspaces, 

X = Hx ®Vx , Y = Hy © Vy , for Vx = Vy = {(j)p) 

and the projection P : Y ^ Hy. Let PFt : Hx —>■ Hy be the projection on Hy of 
the restriction of F to the affine subspace Hx + tfip, PFt{w) = PF{w + tfip). In 
the same fashion, the nonlinearity N -.Y ^ Y gives rise to maps PNt : Hy Hy, 
which we require to be Lipschitz with constant n independent of t G R so that 

[-n,n] n cr(L) = {Ap} . (H) 

The standard Ambrosetti-Prodi map fits these hypotheses. In this case, X cY 
are Sobolev spaces and the derivative /' : R —>• R is bounded by a and b. Set 

7 = (a + b)/2, L = —A — 7 , N{u) = f{u) — ju 

and Ap = Ai, the smallest eigenvalue of —A. Then the Lipschitz constant n of the 
maps PNt satishes n <7 — a = 6 — 7 <A 2 — 7 , so that Ai — 7 < n. 

Theorem 3 Let F : X ^ Y satisfy (H) above. Then for each f G R, the map 
PFt is a bi-Lipschitz homeomorphism, and a diffeomorphism if F is . The 

Lipschitz constants for PFt and {PFt)~^ are independent oft. 

Proof: The proof follows Theorem [1] once the potentially nasty eigenvalue Ap is 
ruled out. Let c be the absolute value of the point in cr(L)\{Ap} closest to 0, so that 
0 < n < c. The operator L ■. X ^Y restricts to L : Hx —>■ Hy, which is invertible 
self-adjoint, and again L~^ : Hy Hy with |1L“^|| < 1/c. The solutions w G Hx 
of PFt{w) = 5 G Hy solve PLu — PN{u) = Lw — PNt{w) = g tor u = w tfip. 
The solutions w correspond to the fixed points of Cg : Hy — >■ Hy, where 

Cg{z) = PNt{L~^z) + g, for Lw = z G Hy. 

The map Cg is a contraction with constant bounded by n/c < 1 (independent of t). 
Now follow the proof of Theorem |TJ ■ 
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The attentive reader may have noticed that the effect of the nonlinearity N 
along the vertical direction is irrelevant for the construction of fibers. 

The same construction applies when the interval [—n, n] defined by the Lipschitz 
constant n of PNt : Hx Hy interacts with an isolated subset I of cr{L) — more 
precisely, I = [—n, n] fl cr(T) and there is an open neighborhood 17 of / C M for 
which I = U n (j{L). In this case P is the orthogonal projection on I, which takes 
into account possible multiplicities. In the special situation when I consists of a 
finite number of eigenvalues (accounting multiplicity), we refer to finite spectral 
interaction between L and N. 

We concentrate on the case when I = {Ap} consists of a simple eigenvalue. 
A more careful inspection of the constants in the Banach iteration in the proof 
above yields the following result ([5], [M]). The image under F of horizontal affine 
subspaces of X are sheets. The inverse under F of vertical lines of Y are fibers. 

Proposition 1 If F is , sheets are graphs of maps from Hy to {(pp) and 
fibers are graphs of maps from {fip) to Hx- Sheets are essentially flat, fibers are 
essentially steep. 

We define what we mean by essential flatness and steepness. Let v{y) be the 
normal at a point y S T of (the tangent space of) a sheet, and t(u) be the tangent 
vector at M S A of a fiber. Then there is a constant e G (0,7r/2) such that pp makes 
an angle less than e (or greater than it — e, due to orientation) with both vectors. 

4 Adapted coordinates and a plan 

Suppose L and N interact at a simple eigenvalue Xp. Write 

F{u) = PF{u) + {F{u), 4>p)4>p = PF{u) + h{u)(j)p 

where the map h : A —>■ M is called the height function. In the diagram below, 
invertible maps are bi-Lipschitz ([34]) or (7^ diffeomorphisms, depending if PFt is 
Lipschitz or . The smoothness of h and = h o $ follow accordingly. 

A = Hx © Vx - 1 Y = Hy © Vy 

<S>-^ = {PFtJd) \ Z' F“=Fo$=(/d,/i“) 

A 

The map F has been put in adapted coordinates by the change of variables 4>: 

{z,t)^{z,h'^{z,t)) . 

Notice that fibers of F are taken to vertical lines in the domain of = A o $. 
Explicitly, the vertical lines {(zo,t) : t G R} parameterized by zq G Hy correspond 
to fibers u{zo,t) = {PFt)~^{zo) + tfip = w{zo,t) + tfip. Thus is just a rank one 
nonlinear perturbation: 

F°'{z, t) = {z, h°‘{z, t)) ^ z + h°‘{z + t(j>p)(j>p . 

In a very strict sense, this is also true of F. In order to make F similar to an 
Ambrosetti-Prodi map, define G = o (—A) : X ^ Y: 

u tA z + t(j)i z + tfii + {h°‘{z + tfii) — t)(j)i = —Au + tp{u)(j)i, 
for some nonlinear functional tjj. We generalize slightly. 
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Proposition 2 Let N he a map. Say L and N interact at a simple eigenvalue 
Xp and L is invertible. Then, after a change of variables, the function 
F = L — N:X^Y becomes G : X ^ Y, G = L + for some : X ^ M.. 

For Ambrosetti-Prodi operators F(u) = —Au — f{u), the nonlinear perturbation 
is given by a Nemitskii map u f(u). It is not surprising that once we enlarge 
the set of nonlinearities new global folds arise. For a map F given in adapted 
coordinates by F°'{z,t) = {z,h°‘{z,t)), appropriate choices of the adapted height 
function h°‘ yields all sorts of behavior. 

The critical set ot F : X ^ Y is compatible with fibers as follows ([J], [?]). 

Proposition 3 Suppose the G^ map F : X ^ Y admits fibers. Then ug is a critical 
point of F if and only if it is a critical point of the height function h along its fiber, 
or equivalently of the adapted height function . 

Isolated local extrema have to alternate between maxima and minima. In par¬ 
ticular, given the appropriate behavior at infinity at each fiber and the fact that all 
critical points are of the same type, we learn from a continuity argument that the 
full critical set G is connected, with a single point on each fiber m)- 

The study of a function F : X ^ Y reduces to three steps: 

1. Stratify X into fibers. 

2. Verify the asymptotic behavior of F along fibers. 

3. Classify the critical points of the restriction of F along fibers. 

The following result is natural from this point of view ([?]). Let F : X ^ Y 
satisfies {H) of Section|31 so that, by Proposition [TJ X stratifies in one dimensional 
fibers {u{z, t) : t £ R}, one for each z e Hy. 

Proposition 4 Suppose that, on each fiber, 

lim {F(u(z,t)),<f>p) = lim h(u{z,t)) = —oo . 

t—^dzoo i—>-±oo 

Suppose also that each critical point of h restricted to each fiber is an isolated local 
maximum. Then F : X ^ Y is a global fold, in the sense that there are homeomor- 
phisms on domain and image that give rise to a diagram as in Theorem\^ 

To verify that such limits exist, one might check hypotheses (F±) in Section 
nn but there are alternatives. Similarly, there are ways of obtaining fibers which 
do not fit the construction presented in Section [3] (this is the case for perturbations 
of non-self-adjoint operators. Section IK31) . The upshot is that there is some loss in 
formulating the three step recipe into a clear cut theorem. 

As trivial examples, h°‘{z,t) = —t^ is a global fold, whereas h°-{z,t) = t^ — t 
has a critical set consisting of two connected components having only (local) folds 
(from Section 1 7. II) . More complicated singularities require the dependence on z: 
not every hber of F (equivalently, vertical line in the domain of F“) has the same 
number of critical points close to a cusp, for example. The reader is invited to check 
that {z,t) !->• {z,t^ — {z,(j))t) is a global cusp, for f any fixed vector in Fly. Higher 
order Morin singularities, considered in Section [71 are obtained in a similar fashion. 
From the Proposition [71 changes of variables on such maps yield nonlinear rank one 
perturbations of the Laplacian which are globally diffeomorphic to the standard 
normal forms of Morin singularities. 
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We consider the standard Ambrosetti-Prodi scenario in the light of this strategy. 
For the function F{u) = — Ait — /(it) defined in Theorem[51 elliptic theory yields all 
sort of benefits — the smallest eigenvalue of the Jacobian DF(u) is always simple, 
the ground state may be taken to be a positive function in X. 

The hypotheses required for the construction of fibers in Theorem|3]do not imply 
the simplicity of the relevant eigenvalue: there are examples for which there is no 
naturally defined functional Ap : X —>• R because two eigenvalues collide. One 
might circumvent this difficulty by forcing the nonlinearity N to be smaller, but it 
turns out that this is not necessary. The hypotheses instead imply the simplicity 
of Ap in an open neighborhood of the critical set C of F, and this is all we need, as 
we shall see in Section |51 

The positivity of the ground state and the convexity of the nonlinearity / are 
used in a combined fashion in the Ambrosetti-Prodi theorem to prove that along 
fibers the height function only has local maxima. Clearly, this is a property only 
of critical points. On the other hand, the nonlinearity N{u) = f{u) is so rigid that 
the standard hypothesis of convexity of / is essentially necessary, as shown in [?]. 
More general nonlinearities require a better understanding of the singularities. 

We now provide more technical details on each of the three steps. 


5 Obtaining fibers in other contexts 

For starters, what if L is not self-adjoint, or X is not Hilbert? 

5.1 Podolak’s approach 

Suppose momentarily that X and Y are Banach spaces. Let L : X ^ Y he a. 
Fredholm operator of index zero with kernel generated by a vector (px and let py 
be a vector not in RanL. Podolak ([IJ) considered the following scenario, for which 
she obtained a lower bound on the number of preimages for a region of Y of vectors 
with very negative component along py. Split X = Fly © Vx where Vx = {<Px) 
and Hx is any complement. Also, split Y = Hy © Vy where Hy = RanL and 
Vy = {(j)y). In particular L : Fix —>■ Hy is an isomorphism. Also, define the 
associated projection P : Y ^ Hy. Write u = w + tpx, U = 9 + spy for w € Hx ■ 
The equation F{u) — Lu — N{u) = y becomes 

L{w + tpx) — N{w + tpx) = Lw — N{w + tpx) = <7 + spy, 

and, as in Theorem [3l we are reduced to solving the map 

Cg : Hy Hy, Cg{z) = PNt{L~^z) + g , for Lw = z G Hy. 

Her hypotheses imply that such maps are contractions. 

5.2 Transplanting fibers 

The estimates arising from spectral theorem in the Hilbert context are easy to obtain 
and possibly more effective. Podolak’s hypotheses are harder to verify. There is 
a possibility: getting fibers in Hilbert spaces and transplanting them to Banach 
spaces. This happens for example when moving from the Ambrosetti-Prodi example 
as a map between Sobolev spaces ([4]) to a map between Holder spaces ([I]). The 
classification of singularities is simpler with additional smoothness (Section [7]). 

Proposition 5 Let F = L — N : X Y satisfy hypothesis {H) of Section [H 
Consider the densely included Banach spaces A C X and B <Z Y allowing for the 
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restriction F : A ^ B for which Vx = Vy A. Suppose that DF{a) : A ^ B 
is a Fredholm operator of index zero for each a £ A. Then fibers of F : X ^ Y 
either belong to A or do not intersect A. 

Said differently, if a point u € X belongs to A then the whole fiber does. 

In the Ambrosetti-Prodi scenario, this proposition seems to be a consequence of 
elliptic regularity, which may be used to prove it. Regularity of eigenfunctions is 
irrelevant: fibers are the orbits of the vector field of their tangent vectors, which are 
inverses of the vertical vector under DF{u), and necessarily lie in A ([?]). Tangent 
vectors are indeed eigenfunctions (j)p(u) of DF[u) at critical points u. 

The fact that sheets and fibers are uniformly flat and steep (Proposition[T]) allows 
one to modify vertical spaces ever slightly and still obtain space decompositions for 
which the Lyapunov-Schmidt decomposition, and hence the construction of fibers 
in Theorem [31 apply. In particular, transplants may be performed even when the 
eigenvector cfp originally used to define the vertical spaces Vx = Vy do not have 
regularity, i.e., do not belong to A G X. We only have to require that A is dense 
in A, so that ftp can be well approximated by a new vertical direction. 

5.3 Fibers and Numerics 

Finite spectral interaction is a very convenient context for numerics. Any question 
related to solving F(u) = g for some fixed g GY reduces to a finite dimensional 
problem in situations of finite spectral interaction, irrespective of additional hy¬ 
potheses. If the interaction involves a simple eigenvalue Ap, one simply has to look 
at the restriction of F to the (one dimensional) fiber associated to the affine vertical 
line through g. 

Smiley and Chun realized the implications of this fact for numerical analysis 
([31], [32]). An implementation for functions F{u) = —Au — f{u) defined on rect¬ 
angles n C was presented in [Sj. In the forecoming sections, we will require 
more stringent hypotheses with the scope of obtaining very well behaved functions 
F — we will mostly be interested in global folds. Such additional restrictions might 
improve on computations, but so far this has not seen to lead to substantial im¬ 
provements on the available algorithms. 




We present an example obtained from programs by Jose Cal Neto (i) and 
Otavio Kaminski. For it = [0,1] x [0, 2], Ai 12.337 and A 2 19.739. Consider 

-Uxx - Uyy - f{u) = g , {x,y) Git , u = 0 in cAI , 

f'(x) = ^ (arctan(^) -‘^x -h Ai , 


/(O)^ 47.12 











9{x, y) = -100(a;(x - l)y‘^{y - 2)) - 35 sin(7ra;) sin(^) . 

On the left, we show the graphs of /', which interacts only with Ai. On the right, 
the height function h associated to the fiber obtained by inverting the vertical fine 
through g. The height value —12.3 is reached by four preimages, displayed below. 
Notice the cameo appearance of the maximum principle: the four graphs sit one on 
top of the other as one goes up along the fiber (this is very specific of interactions 
with Ai of the Laplacian with Dirichlet conditions). 



6 Asymptotics of F on fibers and vertical lines 

We stick to one dimensional fibers and consider two issues. 

1. How does F behave at infinity along fibers? 

2. How do fibers look like at infinity ? 

The first question, to say the very least, is tantamount to characterizing the 
image of F. The second is not relevant for the theoretical study of the global 
geometry of F, since a (global) coordinate system leading to a normal form (like 
{z,t) n- {z,—F) is insensitive to the shape of fibers. On the other hand, for 
numerical purposes, a uniform behavior at infinity of the fibers is informative. 

6.1 F along fibers 

The inverse of a vertical fine zq + Vy, zq G Hy is the fiber u{zo, t) = w(zo, t) + 

F{u{zo, t)) = zo + h‘"(zo, t) (j)p . (*) 
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For a fixed zq G ffy, the map 1 1 —>■ h°'(zo, t) is the adapted height funetion of the 
fiber associated to zq. Clearly, 

/i“(zo, t) = {F{u(zo, t)),4>p) = {L(w{zo, t) + t(j)p) - N{u{zo, t)),4>p) 

so that 

/i“(zo, t) = Xpt - {N{u{zo, t)),(j)p). 

In order to have 


lim {F{u{zo:t)),(j)p) = lim h°‘{zo,t) = -oo 

t^±co i—)-±oo 

and some uniformity convenient to obtain properness as discussed in Section 16.41 
we require an extra hypothesis: 

For each Zq G X, there is a ball U{zq) C X and e, T > 0, c± such that, for 2 G C/(zq), 

{N{u{z,t)),(j)p) > {Xp + e)t + c+ , for f > T, (I/+) 

{N{u{z,t)),(j)p) > {Xp-e)t +c- , ior t <-T. (^-) 

Notice that the asymptotic behavior on each fiber is the same. 


6.2 Asymptotic geometry of fibers 

Again, parameterize fibers as u{z,t) = w{z,t) + t(l)p. Under mild hypotheses, the 
vectors w{z,t)/t have a limit for t ± 00 , which is independent of z. A version of 
this result was originally obtained by Podolak (ED)- 


Proposition 6 Suppose that F : X^Y,F = L — N satisfies hypothesis {H) of 
Section\M Suppose also that, for every u G X, 


PN(tu) 

lim - 

t —^ + 00 't 


NM e u. 


Then there exist w+,W- G Hx sueh that, for every fiber u(z,t) 


V M 

lim- w+\\x = U 

t^+oo " t 


r „w{z,t) 
lim - W- 

f —CX) t 


whieh are respectively the unique solutions of the equations 


= w{z,t) Ftfp, 

\\x = 0 


Lw — PNoc{w + fp) = 0 , Lw + PNao{-w - (j)p) = 0. 

It turns out that Noo = PNoo satisfies the same Lipschitz bound that the func¬ 
tions PNt in Theorem [3J which is why both equations are (uniquely) solvable. 

Fibers are asymptotically vertical if and only if lim|j|_^oo ic(z, t)/t = 0, or equiv¬ 
alently, PNoci^fp) = 0. Indeed, in this case, rc = 0 is the unique solution of 
both equations. This is what happens in the Ambrosetti-Prodi scenario, where 
PNoo{u) = {b — j)Pu~^ — {a — 'y)Pu~ (recall m = — u~), since fp = > 0. 


6.3 Comparing F on fibers and on vertical lines 

One might wish to relate the heights of F along fibers and vertical lines, which are 
easier to handle. In m Podolak presented a scenario in which this is possible. We 
state a version of her result for the case t —>• -l-oo. 

Theorem 4 Let X C Y be Hilbert spaces with X dense in Y. Let L : X ^ Y 
be a self-adjoint operator with 0 G cr(T), a simple, isolated eigenvalue, associated 
to the normalized kernel veetor fp. Set Hy = {fp)'^. Take N : Y ^ Y and 
F = L — N:X^Y so that 
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1. ||iV(w) - 7V(uo)||y < e||u - uoIIy , limj^+oo iV(to)/f = iVoo(u) 

2. {Nooi(l>p),(l)p) = -limt^+oo{F{t(j)p),(j)p)/t > 0 

5. e \\{L\h^)-^\\ < 1/2 , £2 \\{L\Hry"\\ < 1/2 (Noo^p), cj^p) . 


Then, for each fiber {zo,t) in adapted coordinates, 



The number {Nac{(t>p),4>p) gives the asymptotic behaviour of the height of F 
along the vertical line through the origin. The theorem implies that F along the 
upper part of each fiber converges to the same infinity that F along {tcjjp, t > 0}. 

A context in which these hypotheses apply is the Ambrosetti-Prodi operator 
with a piecewise nonlinearity f(u) = (Ap + c)u'^ — (Ap — c)u~ for a sufficiently small 
number c > 0. However, for pairs (Ap — ci, Ap + C 2 ),p ^ 1 in the Fucik spectrum of 
the (Dirichlet) negative second derivative, for which necessarily ci c^ (near Ap), 
the condition involving does not hold and indeed the thesis is not true. 

6.4 Fibers and the properness of F 

From a more theoretical point of view, fibers circumvent the fundamental issue of 
deciding if F is proper. For example ([IS]), the map 


F : (§^) ^ (17° (S^), u^vf F arctan(M) 

is a diffeomorphism from the domain to the open region between two parallel planes, 



Indeed, fibers in this case are simply lines parallel to the vertical line of constant 
functions, and each is taken to such region. 

Perhaps, it would be more appropriate to think of fibers as a tool to show 
properness ([?]). As far as we know, for the Ambrosetti-Prodi map F : Af —>■ F in 
unbounded domains, the properness has been proved only by making use of fibers 
(see Section 0). 

Proposition 7 The map F : X ^ Y satisfying hypotheses (FI) of Section 0 and 
(F±) above is proper if and only if the restriction of F to each fiber is proper. 

Points in the Fucik spectrum of the (Dirichlet) second derivative give rise to 
maps F which take the half-fiber {u{0,t),t > 0} to a single point 0 ([34]), which 
shows that F is not proper, although the image of every vertical line has its vertical 
component taken to infinity. 

A possible definition of a topological degree for F becomes innocuous — the 
relevant information is essentially the asymptotic behavior of F along each fiber. 

7 Singularities 

Generic singularities both of F and of each height function are very special — they 
are Morin singularities. Morin classified generic singularities of functions from R” 
to R" whose derivative at the singularity has one dimensional kernel ([IS]). This 
is sufficient for the study of critical points of height functions on one dimensional 
fibers, by Proposition [3| In order to do the same for the critical points of the 
whole function F : X ^ Y, we need an equivalent classification for singularities of 
functions between infinite-dimensional spaces, which is very similar m, m, ED) 
— this is how we proceed next. 
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7.1 Morin theory in adapted coordinates 

The first step in Morin’s proof makes use of the implicit function theorem to write 
such a singularity at a point (zo,to) in adapted coordinates, as in Section 2) 

: r = iIy e Vy ^ r = ffy e Vy, (z,t) ^ (z,h‘^(z,t)). 

Say is The point (zo,to) is a Morin singularity of order k if and only if 

1 . Dth^izo,to) = ■■■ = D^h‘^{zo,to) = 0 , D'^+^h‘^{zo,to) ^ 0 . 

2. The Jacobian D{h°', Dth°‘,..., D^~^h°‘){zQ, to) has maximum rank. 

Then, in a neighborhood of (zq, to) there is an additional change of variables which 
converts F°‘ to the normal form 

{z,x,f) {z,x,t^^^ +Xit^~^ H- \-Xk-lt) . 

Here the coordinates (z, x) correspond to an appropriate splitting oiY = Y . 

Morin singularities of order 1, 2, 3 and 4 are called, respectively, folds, cusps, 
swallowtails and butterflies. 

Thus, the classification of critical points of F boils down to the study of a 
family of one dimensional maps, the height functions restricted on fibers. The 
first requirement is specific to each fiber (i.e., one checks it for every fixed z near 
zo), whereas the second relates nearby fibers, i.e., one has to change z. Folds are 
structurally simpler than deeper singularities: the behavior along fibers near a fold 
point is always the same — essentially like 1 1—>■ —whereas this is not the case for 
cusps, where close to 1t^ one finds 1t^ ± et. 

There is something unsatisfying in the fact that the relevant properties of the 
critical points of F requires knowledge of some version of the height function. This 
is circumvented by the next result ([?]). 

Proposition 8 Suppose F : X ^ Y is and admits one dimensional fibers. 

Then there is an open neighborhood U of the critical set C with the properties below. 

1. There is a unique map Xp : U ^ M. for which Xp = 0 on C and is an 
eigenvalue of DF elsewhere. 

2. There is a strictly positive function p : U ^ such that 

Xp{u{z,t)) =p{u{z,t)) Dth{u{z,t)) , u(z,t) € U . 

A point uo = u{zo,to) is a Morin singularity of order k of F if and only if 

1. Xp{uo) = • • • = D’f ^Xp{uo) = 0 , D^Xpipio) 0 , 

2. The image of D{Xp ,..., D^~^Xp){uo) together with DtXp{uo) span R". 

There is an analogous characterization in adapted coordinates. 

7.2 Critical points of the height function 

Consider a critical point uo € C C X and the fiber u(zo, t) through it, u(zo, to) = uo. 
From Proposition |S1 uo is a (topological) fold of the height function h restricted to 
the fiber if and only if uo is a topologically simple root of Xp{u) along the fiber, i.e., 
Xp is strictly negative on one side of uo and strictly positive on the other. 
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Once we reduce the issue to checking an eigenvalue along a fiber, derivatives are 
irrelevant: just study the quadratic form of the Jacobian. Clearly, this only handles 
topological equivalence between the function and a fold. 

More explicitly, in standard Ambrosetti-Prodi contexts, Ai(uo) is the minimum 
value of the quadratic form {DF{uo)v,v). The derivative Dtu{zo,to) of the (C^) 
fiber is the eigenfunction (f>i{uo) > 0 , and it is easy to check that Ai increases with 
t by the convexity of the nonlinearity /. This should be compared with differentia¬ 
bility arguments, which require some estimate on (^i(uo) (say, boundedness). 

The fact that all critical points are local maxima for height functions on fibers, as 
required in Proposition |4l suggest hypotheses to be checked only on the critical set 
of F. This is not the case in the original Ambrosetti-Prodi theorem: the statement 
of the theorem has the merit that it makes no reference to the critical set at all, 
an object which in principle is hard to identify. The convexity of the nonlinearity 
handles the difficulty and, rather surprisingly, is essentially necessary ([?]). Further 
examples yielding local maximality are somewhat contrived. 


8 Some examples 

8.1 The non-autonomous case 

The geometric formulation F = L — N is not sufficient to accomodate situations 
of the form F{u{x)) = —Au(x) + f(x,u(x)), the so called non- autonomous case. 
Hammerstein ([IB]) had already considered homeomorphisms of that form. A pos¬ 
sibility is requiring that X and Y are function spaces defined on a domain 12, so 
that the variable x makes sense. The formalism above carries over to this scenario 
without surprises. 

More precisely, as usual X and Y are Hilbert spaces, X dense in Y. The linear 
operator L : X CY ^ P is self-adjoint with a simple eigenvalue Ap associated to a 
normalized eigenvector (f)p. Let P :Y ^ Hy = be the orthogonal projection. 

From the nonlinear term N : D,xY —)• P, define as before PNt : Hy Fly, t gR 
and require a Lipschitz estimate, 

\\PNt{x,wi) - PNt{x,wo)\\y < n\\wi - wqIIf, for wo,wi G Hy, 

so that [—n, n] fl cr{L) = {Ap}, which is the same hypothesis (H) in Section[3l This 
obtains fibers for F : X ^ Y as in Theorem [3l which satisfy the same properties as 
those in the autonomous case, in particular, Proposition [TJ 

The hypothesis which obtain appropriate asymptotic behavior of F along fibers 
are the obvious counterparts of (P-k) and (P—) in Section l6Tl For the classification 
of critical points, we simply do not distinguish between the autonomous and non- 
autonomous case: the subject has become a geometric issue. 

8.2 Schrodinger operators on 

As was surely known by Ambrosetti and Prodi (and [2] is an interesting exam¬ 
ple), the Laplacian with Dirichlet conditions might be replaced by more general 
self-adjoint operators. The approach in this text is flexible enough to handle non¬ 
linear perturbations of Schrodinger operators on unbounded domains yielding global 
folds. In our knowledge there are no similar results in the literature. Tehrani ([33]) 
obtained counting results for Schrodinger operators in R." in the spirit of those 
obtained by Podolak ([ 21 ]), indicated in Section O- 

We state the by now natural hypotheses. Here Y = L^(]R"'). 

1. The free operator T = —A + v{x) : A C P —>■ P is self-adjoint, with simple, 
isolated, smallest eigenvalue Ai and positive ground state (/)i. 
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2. F -.X CY ^Y, F{u) =Tu- f{u) is a map. 

3. The function / G C^(R) satisfies /(O) = 0, M > /" > 0, /'(K) = (a,5) and 
a < Ai < 5 < min{(T(r) \ {Ai}}. 

4. The Jacobians DF{u) ■. X ^ Y are self-adjoint operators with eigenpair 
(Ai(u), (/)i(u)) sharing the properties of (Ai,(/)i). 

Theorem 5 Under these hypotheses, the map F : X ^ Y is a global fold. 

Such hypotheses are satisfied for v{x) = x^/2, the one dimensional quantum 
harmonic oscillator, as well as for the hydrogen atom in for which v{x) = — l/|a:|. 

Hypotheses on the potential of a Schrodinger operator in order to obtain such 
properties are commonly studied in mathematical physics. The interested reader 
might consider H, m, m- More about this in [?]. 

8.3 Perturbations of compact operators 

We recall Mandhyan’s second example of a global fold ([22]), or better, a special 
case of the extension given by Church and Timourian (in])- 

For H C M" a compact subset, let X = C^{V,) and define the compact operator 

K-.X^X, K{u){x)= [ k{x,y)u{y)dy 

Jn 

where the kernel k G C^(fl x fl) is symmetric and positive. Let yi > p ,2 be the 
largest eigenvalues of K. Now let / : R —?> R be a strictly convex function 
satisfying 

0 < lim f{x) < l/yi < lim f'{x) < l/\p 2 \ ■ 

X—¥ — 00 X—^OO 


Theorem 6 Under these hypotheses for K and f, the map 

G : X ^ X, G{u){x) = u{x) — Kf{u{y)) 


is a global fold. 

This is the kind of nonlinear map obtained if one started from the Ambrosetti- 
Prodi original operator F{u) = —Au — f{u) and inverted the Laplacian. Actually, 
one could take another track: instead of inverting the linear part, one might consider 
the inversion of the nonlinear map u f(u), since f is bounded away from zero. 
For maps G{u) = Ku — f{u) obtained this way, we handle the case when AT is a 
general compact symmetric operator K. 

More precisely, let H C R", B = G'^{Uf) and Y = Let K : B ^ B and 

K : Y —>■ F be compact operators which preserve the cone of positive functions. 
Also, K : Y ^ Y has simple largest eigenvalue Xp = \\K\\ and second largest 
eigenvalue As . Let / : R —>■ R be a strictly convex function, with /(O) = 0 if H 
is unbounded. Suppose 

As < a = lim f'(t) < Xp < b = lim f'(t) . 

t—¥ — oo t—^oo 

Theorem 7 The map F : B ^ B , F{u) = Ku — f{u) is a global fold. 

The reader should notice that F is Lipschitz but not differentiable as a map 
from Lf{Tl) to itself. Still, the direct construction of fibers in G^{Vl) is not a simple 
matter, because properness of F is not immediate. Transplanting fibers in this 
example is convenient, and was also used in Mandhyan’s context. 
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8.4 Folds as perturbations of non-self-adjoint operators 

McKean and Scovel m, m) studied the Riccati-like map on functions 

u G L2([0, 1]) ^ w + {D,)-^f{u) G L^iO, 1]), fix) = xy2, 

where (- 02 )”^ is the inverse of the second derivative acting on W^’^i[0,l]) and 
showed that the critical set consists of a countable union of (topological) hyper¬ 
planes. Church and Timourian ([8]) showed that the restriction of such map to 
a neighborhood of one specific critical component is (after global homeomorphic 
change of variables) a fold. The techniques employed are in the spirit of the original 
Ambrosetti-Prodi paper. 

Fibers were relevant in f|25j'). where perturbations of first order differential equa¬ 
tions (clearly, non-self-adjoint operators) were shown to be global folds. An example 
is the map on periodic functions with (generic) convex nonlinearities /, 

F ■.C^iS^)^C°iS^), u^u' + fiu). 

McKean and Scovel ([13]) and Kappeler and Topalov ([13]) considered the same 
map among Sobolev spaces, the celebrated Miura map, used as a change of variables 
between the Korteweg-deVries equation and its so called modified version. 

More recently, a perturbation of a non-self-adjoint elliptic operator (as in [3], 
but with Lipschitz boundary) has been shown to yield a global fold ([30]). 

8.5 Acknowledgements 

The first author thanks the Departamento de Matematica, PUC-Rio, for its warm 
hospitality. The second and third authors gratefully acknowledge support from 
CAPES, CNPq and FAPERJ. We thank Jose Cal Neto and Otavio Kaminski for 
the numerical examples. 


References 

[1] A. Ambrosetti and G. Prodi, On the inversion of some differentiable mappings 
with singularities between Banach spaces, Ann. Mat. Pura Appl. 93 (1972) 
231-246. 

[2] H. Berestycki, Le nombre des solutions de certain problemes semi-lineaires 
elliptiques, J. Eunct. Anal. 40 (1981), 1-29. 

[3] H. Berestycki, L. Nirenberg and S.R.S. Varadhan, The principal eigenvalue 
and maximum principle for second-order elliptic operators in general domains. 
Comm. Pure Appl. Math. 47 (1994) 47-92. 

[4] M.S. Berger and E. Podolak, On the solutions of a nonlinear Dirichlet prob¬ 
lem, Indiana Univ. Math. J. 24 (1974) 837-846. 

[5] M.S. Berger and P.T. Church, Complete integrability and perturbation of a 
nonlinear Dirichlet problem. L, Indiana Univ. Math. J. 28 (1979), 935-952. 

[6] M.S. Berger, P.T. Church and J.G. Timourian, Eolds and cusps in Banach 
spaces, with applications to nonlinear partial differential equations. L, Indiana 
Univ. Math. J. 34 (1985), 1-19. 


15 


[7] F.A. Berezin and M.A. Shubin, The Schrodinger Equation (1991), Kluwer, 
Dordrecht. 

[8] J.T. Cal Neto and C. Tomei, Numerical analysis of semilinear elliptic equa¬ 
tions with finite spectral interaction, J.Math.Anal.Appl. 395 (2012) 63-77. 

[9] M. Calanchi, C. Tomei and A. Zaccur, Fibers and global folds in infinite 
dimension, in preparation. 

[10] P.T. Church and J.G. Timourian, Global fold maps in differential and integral 
equations, Nonlinear Anal. 18 (1992) , 743-758. 

[11] P.T. Church, E.N. Dancer and J.G. Timourian, The structure of a nonlinear 
elliptic operator, Trans. Amer. Math. Soc. 338 (1993) 1-42. 

[12] P.T. Ghurch and J.G. Timourian, Global Structure for Nonlinear Operators in 
Differential and Integral Equations I. Folds, Topological Nonlinear Analysis, 
II, Progr. Nonlinear Differential Equations Appl., 27, Birkhauser, Boston, 
MA, (1997) 109-160. 

[13] P.T. Ghurch and J.G. Timourian, Global Structure for Nonlinear Operators in 
Differential and Integral Equations II.Cusps, Topological Nonlinear Analysis, 
II, Progr. Nonlinear Differential Equations Appl., 27, Birkhauser, Boston, 
MA, (1997) 161-246. 

[14] D.G. Costa, D.G. Figueiredo and P.N. Srikanth, The exact number of solu¬ 
tions for a class of ordinary differential equations through Morse index com¬ 
putation, J. Diff. Eqns. 96 (1992) 185-199. 

[15] C.L. Dolph, Nonlinear integral equations of the Hammerstein type, Trans. 
AMS 66 (1949) 289-307. 

[16] A. Hammerstein, Nichtlineare Integralgleichungen nebst Anwendungen, Acta 
Math. 54 (1929) 117-176. 

[17] T. Kato, Perturbation Theory for Linear Operators, Classics in Mathematics 
(1980), Springer, New York. 

[18] T. Kappeler and P. Topalov, Global fold structure of the Miura map on T^(§), 
Int. Math. Res. Not. 39 (2004), 2039-2068. 

[19] A.C. Lazer and P.J. McKenna, On the number of solutions of a nonlinear 
Dirichlet problem, J. Math Anal. Appl. 84 (1981) 282-294. 

[20] E.H. Lieb and M. Loss, Analysis, Graduate Studies in Mathematics 14 (2001), 
AMS, Providence. 

[21] 1. Mandhyan, Examples of global normal forms for some simple nonlinear 
integral operators. Nonlinear Anal. 13 (1989),1057-1066. 

[22] 1. Mandhyan, The diagonalization and computation of some nonlinear inte¬ 
gral operators. Nonlinear Anal. 23 (1994), 447-466. 

[23] H.P. McKean and J.G. Scovel, Geometry of some simple nonlinear differential 
operators, Ann. Scuola Norm. Sup. Pisa Cl. Sci. 13 (1986), 299-346. 

[24] A. Manes and A.M. Micheletti, Un’estensione della teoria variazionale classica 
degli autovalori per operator! ellittici del secondo ordine, Boll. Unione Mat. 
Ital.. 7 (1973) 285-301. 


16 



[25] I. Malta, N. Saldanha and C. Tomei, Morin singularities and global geometry 
in a class of ordinary differential operators, Topol. Meth. Nonlin. Anal. 10, 
137-169, 1997. 

[26] B. Morin, Formes canoniques des singularites d’une application differentiable, 
C. R. Acad. Sci. 260 (1965), 5662-5665, 6503-6506. 

[27] E. Podolak, On the range of operator equations with an asymptotically non¬ 
linear term, Indiana Univ. Math. J. 25, 12 (1976) 1127-1137. 

[28] B. Ruf, Singularity theory and bifurcation phenomena in differential equa¬ 
tions, Topological Nonlinear Analysis, II, Progr. Nonlinear Diff. Eqs. AppL, 
27, Birkhauser, Boston, MA, (1997) 315-395. 

[29] M. Reed and B. Simon, Fourier Analysis, Self-Adiointness,Academic Press, 
New York (1975). 

[30] B. Sirakov, C. Tomei and A. Zaccur, Global folds given by perturbations of 
non-self-adjoint elliptic operators, in preparation. 

[31] M.W. Smiley, A finite element method for computing the bifurcation function 
for semilinear elliptic BVPs, J. Comput. Appl. Math., 70 (1996) 311- 327. 

[32] M.W. Smiley and C. Chun, Approximation of the bifurcation equation func¬ 
tion for elliptic boundary value problems, Numer. Meth. Partial Diff. Eqs. 16 
(2000) 194-213. 

[33] H.T. Tehrani, A multiplicity result for the jumping nonlinearity problem, J. 
Diff. Eqs. 188 (2003) 272-305. 

[34] C. Tomei and A. Zaccur, Geometric aspects of Ambrosetti-Prodi operators 
with Lipschitz nonlinearities. Analysis and topology in nonlinear differential 
equations, Progr. Nonlinear Differential Equations Appl., 85, Springer, Cham, 
(2014) 445-456. 


Marta Calanchi, Dipartimento di Matematica, Universita di Milano, 

Via Saldini 50, 20133 Milano, Italia 

Carlos Tomei and Andre Zaccur, Departamento de Matematica, PUC-Rio, 
R. Mq. de S. Vicente 225, Rio de Janeiro, RJ 22453-900, Brazil 

mart a. calanchi® unimi. it 
carlos.tomei@gmail.com 
zaccur. andre@gmail. com 


17 



